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Abstract — This  paper  presents  a  zone  model  of  liquid  pool  tires  burning  in  an  open  environment. 
The  model  is  a  system  of  coupled  non-linear  equations  which  contains  physical  parameters  for 
which  reliable  values  are  not  currently  available.  With  a  judicious  estimate  of  these  parameters, 
the  model  was  tested  on  arious  size  fires  burning  three  types  of  fuel.  The  predicted  results  are 
comparable  with  experimental  results  from  three  groups  of  earlier  workers.  Sensitivity  analyses 
were  run  to  determine  the  relative  influence  of  each  estimated  parameter  on  the  predicted  results. 
These  analyses  show  that  efficiency  of  combustion  in  the  fire  is  the  most  important  estimated 
parameter  in  this  model. 


INTRODUCTION 


With  the  rising  problem  of  plane  crash  fires  and  accidental  fuel  spill  fires  on  ships  decks, 
there  is  need  to  understand  the  burning  characteristics  of  these  fires  in  order  to 
develop  effective  ways  of  fighting  them.  This  paper  presents  a  simple  technique  for 
predicting  these  characteristics  in  freely  burning  pool  fires.  The  majority  of  existing 
fire  models  (e.g..  Quintiere.  1977;  Mitler  and  Emmons,  1980;  Bullen  and  Thomas. 
1979;  Satoh  et  a/.,  1980)  were  developed  principally  to  predict  these  characteristics 
in  confined  fires.  This  is  because  of  the  need  to  establish  more  effective  ways  of 
extinguishing  building  fires  which  occur  more  often  and  hence  cause  most  of  the 
losses  in  property  and  human  lives. 

Theoretical  fire  models  are  often  characterized  as  “field"  or  “zone"  models.  Field 
modeling  yields  point  by  point  solution  of  partial  differential  equations  that  describe 
the  fire  physics.  In  contrast,  the  zone  model  divides  the  system  into  sections  and 
solves  a  set  of  equations  which  describes  what  happens  within  each  section  and  also 
the  interaction  between  sections.  Hence,  in  general,  zone  modeling  yields  less  detailed 
and  less  accurate  results,  although  it  is  an  easier  technique  for  obtaining  practical 
solutions.  Furthermore,  the  accuracy  of  any  such  solution  depends  on  (a)  how  well 
the  equations  describe  what  happens  within  and  between  zones,  and  (b)  the  value  of 
physical  parameters  (e.g.,  fuel  surface  emissivity)  contained  in  the  equations. 
Measured  values  for  some  of  these  parameters  are  currently  not  available. 

The  objectives  of  this  study  are  (I)  to  develop  a  zone  model  for  predicting  the 
steady  state  characteristics  of  various  liquid  pool  fires  in  the  open,  and  (2)  to  determine 
through  a  sensitivity  analysis  the  relative  effect  of  any  error  in  the  estimated  values  of 
the  required  physical  parameters.  The  results  of  this  analysis  will  then  serve  as  a  guide 
for  adapting  this  model  for  predicting  transient  characteristics  of  other  fires. 

The  next  section  presents  the  model  which  comprises  a  system  of  non-linear 
coupled  equations  which  is  solved  by  a  hybrid  numerical  algorithm  developed  by 
Powell  (1970).  With  an  estimated  set  of  physical  parameters,  the  model  was  tested  on 
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various  size  pool  fires.  Finally,  a  sensitivity  analysis  was  undertaken  to  determine 
how  an  error  in  the  estimate  of  each  parameter  would  affect  the  predictions  of  the 
model. 


FIGURE  1  Schematic  of  a  liquid  pool  fire  showing  the  various  zones  and  mass  flows. 


MATHEMATICAL  MODEL 

The  model  considers  a  liquid  pool  fire  burning  freely  in  a  wind-free  open  environment. 
The  system  is  divided  into  three  zones;  namely,  fuel,  combustion  zone  (fire)  and 
plume  (see  Figure  I).  Each  zone  is  assumed  to  have  uniform  properties  and  the  model 
accounts  for  heat  and  mass  transfer  between  zones. 

The  fuel  is  contained  in  a  circular  pan  of  diameter  d,  with  the  liquid  level  kept 
constant  and  flush  with  the  top  of  the  pan.  Fuel  is  assumed  to  evaporate  from  the 
surface  at  the  liquid  boiling  point,  Ts,  while  the  bulk  of  the  liquid  remains  at  Tb ■  Hence 
the  fuel  evaporation  rate  is 

[Ashs(Tf  —  Tk)  +  ilrs.net] 

mv  = -  (1; 

A/f„  +  Cfurl  (Tk  -  Tb) 

where  qm.ntt— the  net  radiant  heat  gain.f  A  heat  conduction  term  (Hottel,  1959) 
has  been  neglected  in  Eq.  (I).  This  introduces  little  error  for  the  pan  diameters 
(  •  15  cm)  of  interest  here. 

As  the  fuel  vapor  rises,  it  entrains  the  surrounding  air  and  burns.  Part  of  the  heat 
generated  is  transferred  back  to  the  fuel  surface  where  it  is  used  to  evaporate  more 
fuel  and  thus  keep  the  fire  burning.  The  resulting  flame  is  assumed  to  be  turbulent 

tThe  symbols  are  explained  in  the  Nomenclature. 
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with  negligible  entrainment  in  the  vertical  direction.  For  simplicity  the  flame  is 
divided  into  two  zones;  namely,  (a)  the  combustion  zone  where  all  the  combustion 
reactions  are  assumed  to  take  place,  and  <b)  the  plume  which  is  comprised  of  hot 
flue  gases.  Therefore,  the  heat  balance  equation  in  the  combustion  zone  is: 


y<jen .  '  Maf  C’aTa  M  v  C r  7*  - 

nhxCgTf  AJiATf  -  T*)  ■  qrf.net 

(2) 

where  the  total  mass  flow  out  of  the  combustion  zone  is 

ii  is* 

m«r  ■  m, : 

(3) 

the  heat  generated  in  the  fire  is 

q<jen. 

j 8mv\Hr 

(4) 

and  the  net  radiant  heat  loss  from  the  fire.  Note  that  a  parameter  fi  has  been 

introduced  to  account  for  combustion  efficiency,  i.e..  the  fraction  of  the  heat  of 
combustion  that  is  actually  released  in  the  fire. 

The  rate  of  air  entrainment  into  the  combustion  zone  is  expressed  in  terms  of 
height  above  the  fuel  surface,  Hf,  Froude  number  Fr ,  heat  generation  in  the  fire  and 
properties  of  the  ambient  atmosphere  (Fang,  1973); 


where 


M  I/  A//c.voA 
Mali  :  yCaTtt  ) 


and  the  normalized  density  at  the  fuel  surface  is 


Po  M/uei  7-00/(29  T„). 


(5) 


(6) 


(7) 


Fang  (1973)  has  suggested  that  the  effects  of  radiation  in  Eq.  (6)  be  approximated  by 
reducing  the  value  of  SHr-  Based  on  the  values  of  radiative  heat  fraction  reported 
in  the  literature  (e.g..  Modak  and  Croce.  1977),  we  chose  to  multiply  AA/<-  by  a  factor 
of  0.7.  Note  that  Eq.  (7)  implies  that  the  mass  fraction  of  fuel  vapor  at  the  fuel  surface 
is  I.  Finally,  the  remaining  terms  in  Eq.  (5)  are 


and 


where 


Hf  .»■„  1.49  0,'  »Z*5; 


,(4 .  J-V 

V  P  0  .Vo,  ' 

a/*(\  —  oj) 


(8) 

(9) 


(10) 
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2.16 

and  the  modified  Froude  number, 

Z  =  p,,' />' *.  (1!) 

Any  unreacted  fuel  and/or  oxygen  leaving  the  combustion  zone  mixes  with  com¬ 
bustion  products  and  as  this  hot  mixture  rises  it  entrains  more  of  the  surrounding  air 
which  cools  it  down.  Thus,  the  energy  equation  in  the  plume  zone  is 

lilxsCffT/  ■  WapCaTa  --  thpCpTp  •  qrp.net*  (12) 

where  the  total  mass  leaving  the  plume  is  given  by 


riit >  *  w»«  f  ritap 


(13) 


and  qrp.ntt^-  the  net  radiant  heat  loss  from  the  plume.  The  air  entrainment  rate  in 
the  plume  region  is  given  by 


Map 


mSs 

Pss' 


4-  I 


-  I 


(14) 


where  the  density,  radius  and  Froude  number  of  the  gas  column  at  the  combustion 
zone/plume  interface  are,  respectively 


Pss 


(15) 


a/Fr  \ 1/5  T(  I  +  y/.vo;)5 


/  a/Fr  \15p 
.v„«  =  .Vo  1.191—— - 

\COl,5(l  —  Co)/ 


—  +  y/xo2 ' 

Po 


1  10 


(po')2,: 


(16) 


and 


F,s  -  (—-)7 

\pss  Pair'  I 


(17) 


Furthermore.  Fang's  analysis  has  shown  that  the  visible  flame  height  correlates 
reasonably  well  with  the  point  on  the  plume  axis  where  the  temperature,  7),  falls  down 
to  about  700:K.  Presently,  the  current  model  uses  Fang’s  expression  for  the  flame 
height,  thus. 


Hp  --  v0(  1  -49  4-  0.916  02' *)  0i>  5Z25,  (18) 


where 
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Tx 

7V  - ;  T,  *  7(XTK. 

T-r 


(19) 


(20) 


However,  other  expressions  for  flame  height  (e.g.,  Zukoski.  1981)  can  equally  be 
used  in  the  model. 

Finally,  the  net  radiant  heat  transfer  between  the  three  zones  is  derived  by  making 
the  following  assumptions: 

a)  The  combustion  zone  and  plume,  each  is  a  grey  absorbing  and  emitting  medium 
with  homogeneous  properties. 

b)  Emissivity.  t,  is  given  by 

e  —  1  —  e  (21) 


where  the  mean  beam  length  is  expressed  as  (Orloff.  1980) 


3.6  V 


(22) 


and  x  is  the  absorption  coefficient. 

c)  The  volumes  of  the  combustion  zone  and  the  plume  are  approximated  by 
frustrums  of  right  cones,  with  dimensions  obtained  from  Fang's  analysis. 

d)  Transmissivity  of  the  fuel  surface  is  zero. 

Thus,  the  radiant  heat  transfer  at  the  fuel  surface,  in  the  combustion  zone  and  in  the 
plume  are  given,  respectively  by  (Eckert  and  Drake,  1972). 


4r*,net  '  " -  [&•  Cft*];  (23) 

I  —  f« 

t jr/.tift  A sFf(Chf  B , )  A/f  f( <’ b<  —  (' hfi)  A  p  Ff{Cbj  t?bp)  (24) 

and 

</rp,ne'  A  ft  Fp(ehp  —  —  BnTfF*p)  •  ( Ap  •  tfY  p^F  p(ehp  —  €  ha)  i  (25) 

wheie 

e*  -  oT *  (26) 


is  the  black  body  emissive  power: 

B*  Fi?bt  (I  ~  f*){r/es/  Tf[FP€hpFnp  esp(  1  —  Ftp)] } 


(27) 
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is  the  fuel  surface  radiosity  and  the  view  factor  between  the  plume  and  surface  is 

128) 


(29) 


[Hf  1 

,  If  Hr  1  ,1 

2  4  2 

Fhp  0.5 

.  ,1V-2 

'-lb.;2  1 

yJ\ 

Furthermore,  the  radius  of  the  top  of  the  plume  is  given  by 
vP  1 2*p(H„  -  Hf)  r«,  I.I031  5Z2\ 

where 


#3 


■  — )’ 
\  po  .Vo,  / 


‘•’H I  —  p*,') 


(30) 


Finally,  the  areas  of  the  plume  surface,  the  lateral  surface  of  the  combustion  zone 
and  the  plume/fire  interface  are  given,  respectively  by 

AP  M.w,  ■  yp)[(H„  -  Hf)*  <y„  -  tv-)4]1  (31) 

A,  -  Myo  •  iv.)  [Hr  v  (.t  o  -  i •«)*!>  2  (32) 

and 

Ap  7r\  ',.<2.  (33) 


NUMERICAL  SOLUTION 

The  system  of  non-linear  algebraic  equations  derived  in  the  previous  section  can  be 
put  in  the  form 

,r2  .  . .  x„)  ~  0,  i  -  1.2...  ii.  (34) 

The  set  was  then  solved  using  a  computer  algorithm  developed  by  Powell  (1970). 
This  algorithm  is  a  variation  of  the  classical  Ncwton-Raphson  iteration  such  that  its 
fast  convergence  is  maintained,  while  the  problems  which  arise  in  the  iteration  when 
the  Jacobian  becomes  singular  or  when  starting  from  a  poor  initial  guess  of  the 
solution  are  avoided. 

The  current  computer  code  was  written  such  that  only  three  equations  [Eqs.  ( 1 ), 
(2)  and  (12)]  in  the  set  are  actually  solved  with  Powell's  algorithm  for  the  three 
unknowns  fhv,  Tf  and  Tv.  However,  at  each  iteration  current  values  of  mr,  Tf  and 
Tp  are  used  to  evaluate  the  rest  of  the  equations.  This  is  equivalent  to  solving  these 
equations  exactly  at  each  iteration.  This  arrangement  is  possible  because  one  can  in 
principle  write  the  above  three  equations  explicitly  in  terms  of  mr,  Tf  and  Tp  using  the 
other  equations.  This  technique  saved  considerable  computer  time. 

In  addition  to  the  initial  guess  of  the  solution  of  the  system  of  equations.  Powell's 
algorithm  requires  the  user  to  supply  other  numerical  parameters;  the  key  ones  being 
the  initial  iteration  step  size  and  a  convergence  criterion.  The  iteration  is  said  to 
converge  when 
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In  the  double  precision  version  of  the  current  code  where  it 


RESULTS  AND  DISCUSSION 

The  model  has  been  tested  on  gasoline,  kerosene  and  methanol  pool  fires  of  various 
sizes.  Pan  diameter  was  varied  between  15  and  200  cm,  a  range  which  covers  various 
stages  of  turbulence  in  the  fire.  The  properties  of  the  three  fuels  used  in  the  calcu¬ 
lations  are  given  in  Table  I.  One  recognizes  that  gasoline  and  kerosene  are  both 
variable  mixtures  of  hydrocarbons  and  hence  their  properties  will  depend  on  the 
mixture.  However,  for  the  purpose  of  this  work,  handbook  (Bolz  and  Tuve,  1976) 
values  of  these  properties  have  been  used. 

Table  II  lists  the  estimated  values  of  the  physical  parameters  for  the  three  liquid 
pool  fires.  These  estimates  are  close  to  the  values  used  in  related  works  (Quintiere, 


TABLE  I 

Properties  of  the  fuels  used 


Property 

Gasoline 

Fuel 

Kerosene 

Methanol 

Specific  heat  (liquid),  C/„t i  W-sec/kg  K 

2364 

2090 

2540 

Boiling  point,  T,  K 

410.8 

499.67 

337.8 

Stoichiometric  air/fuel  ratio 

14.97 

14.99 

6.47 

Heat  of  combustion.  A//<-  W-sec/kg 

4.39  x  107 

4.30  x  107 

2.23  x  107 

Heat  of  vaporization,  -iff,  W-sec/kg 

2.696  x  I05 

2.0  x  10’ 

11  x  I05 

Molecular  weight,  Mr„,i 

126 

154 

32 

Liquid  density  at  boiling  point,  kg/m3 

684.1 

740.2 

786.5 

TABLE  II 

Best  estimate  of  the  various  parameters 


Parameter 

Value  used  for 
gasoline  and  kerosene 
fires 

Value  used  for 
methanol  fire 

0 

0.73 

0.95 

0.42 

0.42 

h.  (W/m*  K) 

2.5 

2.5 

M/Ma 

3.0 

1.0 

K/  (ITT1) 

1.05 

0.85 

«„(m  ') 

0.75 

0.35 

*/ 

0.25 

0.25 

0.15 

0.15 
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1977;  Miller  and  Emmons,  1981;  Markstein,  1978)  and  they  were  tested  by  comparing 
predicted  fire  characteristics  with  measured  values  (Blinov  and  Khudiakov.  1957). 
Note  that  the  values  of  some  of  the  parameters  for  methanol  fire  (Table  l()  are 
different  from  similar  values  for  kerosene  and  gasoline  fires.  In  general,  methanol 
fires  burn  cleaner  than  kerosene  or  gasoline  fires  of  the  same  size.  Consequently.  (3 
for  the  methanol  fire  was  estimated  close  to  1.0  while  absorption  coefficients  were 
estimated  lower  than  the  values  for  gasoline  and  kerosene  fires.  In  addition,  the 
molecular  weight  of  methanol  is  nearly  equal  to  that  of  air,  hence  the  normalized 
average  molecular  weight  in  the  methanol  fire  was  estimated  to  be  1 .0. 

One  would  expect  some  of  these  physical  parameters  to  vary  with  pan  size.  In  fact, 
recent  work  by  Markstein  (1978)  suggests  that  the  average  fire  temperature  would 
vary  inversely  witi.  pan  size,  while  the  absorption  coefficient  would  increase  with  pan 
size.  Also  in  a  recent  paper  by  Tewarson  et  al.  (1980)  it  was  reported  that  for  small 
size  fires  (d=  10  and  30  cm)  the  combustion  efficiency  does  not  vary  appreciably  with 
d.  Although  this  conclusion  is  consistent  with  the  current  work,  one  would  expect 
that  over  a  wider  range  of  pan  size,  fi  would  vary  inversely  with  d.  These  effects 
probably  result  from  inefficiencies  such  as  soot  formation,  etc.  However,  since  it  is 
not  known  exactly  how  these  parameters  vary  with  pan  size,  ,8  and  MjMa  were  allowed 
to  vary  in  this  work  (for  d ^  100  cm)  in  order  to  obtain  results  that  are  comparable  to 
measured  values.  Hence,  for  pan  diameters  greater  than  100  cm.  /3  was  decreased 
to  0.64  while  /i?/7Wa  was  increased  to  4.0.  One  could  equally  increase  «/  for  large 
fires  to  improve  the  results. 

Table  III  compares  the  predicted  characteristics  of  three  30  cm  diameter  liquid 
pool  fires  (methanol,  gasoline  and  kerosene)  with  the  same  characteristics  measured 
by  three  independent  groups,  namely  Blinov  and  Khudiakov  (1957).  Rasbash  et  al. 
(1956)  and  Burgess  et  al.  (1959).  The  fuel  evaporation  velocity,  t'n.  is  the  rate  at  which 
the  liquid  level  would  fall  if  it  was  not  maintained  flush  with  the  top  of  the  pan. 
Table  III  reveals  that  with  the  given  values  of  the  physical  parameters,  the  current 
model  predicts  the  characteristics  of  unconfined  fires  reasonably  well. 

Next,  the  model  was  further  tested  by  applying  it  to  kerosene  and  gasoline  pan 
fires  of  base  diameters  ranging  from  15  to  200  cm.  The  results  are  summarized  in 


TABLE  III 

Comparison  of  measured  and  predicted  characteristics  for  30  cm  base  diameter  liquid 

pool  fires 


Source 

Fuel 

Evaporation 
velocity  t  o 
(mm/min) 

Temperature 
of  combustion 
zone  7/  K 

Visible 
flame  height. 
HP.  (cm) 

This  study 

Methanol 

0.85 

1504 

74 

Rasbash  et  at. 

Methanol 

0.99 

1491 

-- 

Burgess  et  at. 

Methanol 

0.9 

— 

This  study 

Gasoline 

2.6 

1351 

104 

Rasbash  el  al. 

Gasoline 

2.3 

1299 

Blinov  and  Khudiakov 

Gasoline 

2.5 

100 

This  study 

Kerosene 

1.7 

1318 

88 

Rasbash  et  al. 

Kerosene 

1.55 

1263 

-- 

Blinov  and  Khudiakov 

Kerosene 

1.8 

mi 
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Some  predicted  characteristics  of  various  size  fires 


Characteristics 

Gasoline  fires 

Kerosene  fires 

Pan  diameter  (cm) 

15  50  1.10 

IS  50  130 

Fuel  evaporation  rate, 
kg/m2-sec 

0.0(8 

0.042 

0.047 

0.012 

0.029 

0.033 

T),  °K 

1395 

1325 

1165 

1353 

1294 

1135 

r„,  "K 

525 

512 

562 

522 

510 

559 

kW 

10.44 

262.8 

1771.8 

6.4 

178.9 

1223.0 

qn,  ,„t ,  kW 

2.0 

63.5 

431.6 

1.3 

44.6 

312.1 

qr,.„. u  kW 

0.13 

3.9 

31.0 

0.089 

3.1 

25.3 

A,h,(T/-T,),  kW 

0.04 

0.45 

2.5 

0.038 

0.39 

2.1 

Hp,  m 

0.5 

1.8 

3.76 

0.41 

1.55 

3.27 

Hi,  m 

0.18 

0.67 

1.68 

0.15 

0.58 

1.47 

that,  kg/sec 

0.0069 

0.173 

1.38 

0.004 

0.121 

0.978 

map,  kg/sec 

0.034 

0.855 

4.24 

0.021 

0581 

2.9 

Table  IV  and  also  in  a  i>o  versus  d  plot  (symbols)  in  Figure  2.  The  solid  lines  represent 
a  plot  of  Blinov  and  Khudiakov's  data.  From  Figure  2  and  Table  IV  one  could 
deduce  that  for  small  pan  diameters  ( d 50  cm),  the  current  model  does  very  well  in 
predicting  the  characteristics  of  freely  burning  pool  fires.  Unfortunately,  it  does  not 
do  as  well  for  large  diameter  fires  ( d >  100  cm).  That  can  be  attributed  to  the  follow¬ 
ing:  (a)  since  the  physical  parameters  were  “optimized”  for  small  fires,  some  of  them 
become  less  accurate  as  d  increases;  (b)  the  flame  shape  for  large  fires  is  more  like  a 
cylinder  rather  than  like  two  frustrums  (see  Figure  I)  as  was  assumed  in  the  current 
model;  and  (c)  the  increasing  effects  of  soot  particles  (e.g.,  on  k)  have  not  been 
included  in  the  model. 

Table  IV  presents  some  of  the  predicted  characteristics  of  gasoline  and  kerosene 
fires  of  various  sizes.  Some  of  these  data  reveal  trends  similar  to  what  has  been 


FIGURE  2  Comparison  of  predicted  and  measured  fuel  evaporation  rates  for  various  size  fires. 
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observed  by  earlier  investigators.  For  example,  the  mass  evaporation  rale  and  the 
ratio  of  radiative  to  convective  heat  transfer  to  the  fuel  surface  both  increase  with 
fire  size  (Modak  and  Croce,  1977;  Hottel,  1959).  while  the  average  fire  temperature 
varies  inversely  with  fire  size  (Markstein,  1978). 

Table  IV  also  shows  that  the  current  model  predicts  visible  (lame  heights  which  are 
higher  than  values  reported  by  Blinov  and  Khudiakov.  especially  for  large  fires. 
Actually,  the  visible  flame  height  is  a  very  difficult  quantity  to  measure  since  the  flame 
is  never  steady.  Hence,  it  is  difficult  to  say  how  "inaccurate"  the  predicted  mean 
heights  actually  are.  However,  we  expect  that  the  actual  mean  visible  flame  height 
should  lie  between  Hf  and  Hp. 


SENSITIVITY  ANALYSIS 

As  indicated  in  the  previous  section  the  measured  values  of  some  of  the  physical 
parameters  used  in  the  current  model  are  not  yet  available  in  the  literature.  This  is 
because  the  measurement  of  these  parameters  is  not  only  difficult  but  also  places  high 
demands  on  resources.  Consequently,  there  is  need  to  determine  the  relative  influence 
of  these  parameters  on  the  model  predictions  so  as  to  appropriately  allocate  resources 
for  measuring  them. 

A  sensitivity  analysis  was  undertaken  to  determine  the  relative  influence  of  the 
following  parameters  on  the  model  predictions;  combustion  efficiency,  jS;  air  entrain¬ 
ment  coefficients,  a/  and  ap\  fuel  surface  emissivity,  fuel  surface  connective  heat 
transfer  coefficient,  //»;  normalized  mean  molecluar  weight  in  the  fire.  M;  Mn  and 
absorption  coefficients  */  and  kp.  Two  methods  of  sensitivity  analysis  were  used; 
namely,  a  "brute  force"  method  ( e.g .,  Ndubizu  and  Durbetaki,  1978).  and  a  Fourier 
Amplitude  Sensitivity  Test  (FAST).  With  the  first  method,  each  parameter  was 
varied  independently  while  the  rest  of  the  parameters  were  held  constant  at  their  base 
values.  A  base  value  is  the  best  estimate  of  each  parameter  (e.g..  Table  I).  For  each 
variation  of  a  parameter  a  new  set  of  fire  characteristics  was  calculated.  Considering 
the  problem  of  fire  spread,  the  liquid  fuel  evaporation  rate  is  deemed  the  key  prediction 
of  the  current  model.  Hence,  Figure  3  summarizes  the  results  of  the  analysis  in  a  plot 
of  percent  change  in  parameter  versus  percent  change  in  fuel  evaporation  rate. 

Using  Figure  3,  one  can  rank  the  various  parameters  in  their  order  of  importance 
in  influencing  the  predicted  value  of  the  fuel  evaporation  rate.  This  ranking  would  be: 
combustion  efficiency,  fuel  surface  emissivity,  fire  absorption  coefficient,  normalized 
mean  molecular  weight  of  species  in  the  fire,  air  entrainment  coefficient  in  the  com¬ 
bustion  zone,  fuel  surface  convective  heat  transfer  coefficient,  and  finally  absorption 
coefficient  and  air  entrainment  coefficient  in  the  plume.  It  is  interesting  to  note  that 
the  parameters  for  the  plume  have  no  measurable  effect  on  the  predicted  fuel  evapor¬ 
ation  rate.  This  is  expected,  considering  the  low  temperatures  in  the  plume  and  the 
view  factor  between  the  plume  and  the  fuel  surface. 

The  main  drawback  of  the  above  “brute  force"  method  is  the  assumption  that  the 
parameters  are  independent  and,  hence,  that  they  can  be  varied  one  at  a  time.  This 
drawback  is  avoided  in  the  FAST  method  developed  by  Cukier  et  al.  (1973)  (also  see 
Cukier  el  al.  (1975)].  With  this  method,  the  parameters  A ((i-  1,  2  ...  it)  are  varied 
simultaneously  through  the  introduction  of  a  variable,  s.  such  that 


Ai  A,®»,if(. 


(36) 


%  CHANGE  IN  PARAMETER 

FIGURE  3  Sensitivity  of  the  various  parameters  on  fuel  evaporation  rate. 


whprp 

u,  «,<«>  sin^v.  (37) 

A,<0’  is  the  best  estimate  of  A  >  (e.g..  Table  l).  The  frequency  4n  is  a  positive  integer 
and  ui<0'  is  the  uncertainty  on  A,  which  defines  the  range  between  which  A*  is  varied. 
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Thus  the  predicted  characteristics  (<■.,?..  »i.)  are  periodic  functions  of  v  with  period 
2j r.  Each  characteristic  is  then  Fourier  analyzed  and  the  Fourier  coefficient  for  each 
frequency  corresponds  to  the  sensitivity  (to  the  characteristic)  of  the  parameter  to 
which  that  frequency  was  assigned.  Hence. 


.v 

2  V 

N  lil,  (.V; 

A  _ _ 

j  l 


)  sin  i/i;  sj:  i  1.2..  . .  #i. 


w  hert 


2ny'.V.  /  1.2...  A’. 


(38) 


(29) 


The  frequencies  (I.  21.  31.  37.  45.  49)  and  A'  (  182)  used  in  the  current  analysis  are 

obtained  from  the  work  of  Shaibly  and  Shuler  (1973).  Assignment  of  frequencies 
to  the  various  parameters  was  arbitrary.  However,  a  different  permutation  of  the 
frequencies  did  not  change  the  conclusions  drawn  from  this  analysis  (Shaibly  and 
Shuler.  1973).  The  uncertainty,  i/"".  for  each  parameter  was  chosen  as  0.3.  In  this 
way.  all  the  parameters  varied  within  ranges  that  are  physically  possible.  For  example. 
fi  varied  between  0.99  and  0.54  for  a  30  cm  diameter  pan  gasoline  fire. 

Table  V  presents  the  results  of  the  sensitivity  analysis  using  the  FAST  method. 
For  a  given  fire  size,  the  relative  magnitude  of  the  coefficients  depicts  the  sensitivities 
of  the  various  corresponding  parameters  on  riu-  Incidentally,  the  magnitude  of  the 
coefficients  in  the  two  fires  should  not  be  compared;  the  different  magnitudes  result 
from  the  definition  of  <f>  in  Eq.  (38).  Thus.  Table  V  shows  that  j8  is  the  most  sensitive 
while  h,  is  the  least  sensitive  parameter  in  both  the  large  and  small  size  fire.  In  fact, 
for  a  small  size  fire  (</  30  cm).  Table  V  can  be  used  to  rank  the  parameters  in  exactly 
the  same  order  as  was  obtained  with  the  “brute  force"  method.  This  suggests  that  for 
small  size  fires  the  various  parameters  may  be  assumed  independent  to  a  good 
approximation. 


TABLE  V 

C  oefficients  for  the  sensitivity  of  the  various  parameters  on  fuel 
evaporation  rate 


Coefficient 

Coefficient 

Parameter 

d  30  cm 

( 1  1 30  cm 

Combustion  efficiency,  & 

7.3 

212.7 

Fuel  surface  emissivity,  <« 

3.4 

106.4 

Fire  absorption  coefficient,  k, 

1.54 

34.68 

Normalized  molecular  weight,  f'M„M 

-  0.95 

-  36.75 

Air  entrainment  coefficient,  «, 

0.71 

26.3 

Convective  heat  transfer  coefficient,  h. 

0.48 

7.08 

Considering  the  feedback  mechanism  which  controls  the  physics  of  liquid  or  solid 
fuel  fires,  it  is  not  surprising  that  the  parameters  which  govern  the  amount  of  heat 
liberated  in  the  fire.  (/3)  and  the  amount  fed  back  to  the  fuel  (e*)  should  be  the  con¬ 
trolling  parameters  in  this  model.  Recall  from  Table  IV  that  radiation  is  the  principal 
mode  of  heat  transfer  to  the  fuel  surface.  Furthermore,  both  methods  revealed  the 
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inverse  relationship  between  the  normalized  mean  molecular  weight  in  the  fire  with 
the  fuel  evaporation  rate.  For  a  given  fuel  and  pan  size,  increase  in  means 

increased  proportion  of  higher  molecular  weight  species  in  the  fire.  This  leads  to  a 
reduced  7}  and  consequently  reduced  m ,. 

Finally,  the  results  of  the  sensitivity  analyses  suggest  that  experimentalists  ought 
to  give  priority  to  the  measurement  of  the  efficiency  of  combustion  in  the  fire  rather 
than  any  other  unknown  parameter  needed  in  fire  models.  In  fact.  F.mmons  (1978) 
in  a  recent  paper,  has  indicated  the  need  to  measure  this  parameter.  Unfortunately, 
this  may  not  be  an  easy  task  since  one  would  expect  the  efficiency  of  combustion  in 
the  fire  to  depend  on  (among  other  things)  the  properties  of  the  fuel  (Tewarson  cl  al.. 
1980).  the  chemistry  of  combustion  in  the  particular  fire  and  turbulence. 


CONCLUSIONS 

A  zone  model  has  been  developed  to  predict  the  characteristics  of  liquid  pool  fires  in 
the  open.  This  model  includes  several  physical  parameters,  values  of  which  have  not 
yet  been  measured.  With  a  judicious  estimate  of  these  parameters,  the  model  was 
shown  to  predict  results  which  are  comparable  with  measured  values.  Sensitivity 
analyses  were  run  to  determine  the  relative  influence  of  the  various  (estimated) 
parameters  on  the  predicted  results.  These  analyses  show  that  the  combustion 
efficiency  in  the  fire  is  the  most  critical  estimated  parameter  in  the  model. 

Finally,  with  the  current  estimate  of  these  parameters,  the  model  can  now  be 
adapted  to  predict  transient  characteristics  of  pool  fires  in  a  completely  sealed 
enclosure  where  the  effects  of  oxygen  depletion  and  pressure  build-up  are  very  im¬ 
portant. 


NOMENCLATURE 

A  area  (m-) 

B  radiosity 

C  specific  heat  at  constant  pressure  (fF—  seckg  K.) 
d  diameter  of  fuel  container  (m) 

fire  base  diameter 
t’h  black  body  emissive  power 

F  view  factor 

Fr  Froude  number  (Eq.  8) 

Fss  Froude  number  (Eq.  17) 

g  acceleration  due  to  gravity  (m  sec2) 

lu  convective  heat  transfer  coefficient  ( WJ m2  K.) 

A H,  latent  heat  of  vaporization  of  the  fuel  (fF-sec/kg) 
A  He  heat  of  combustion  of  the  fuel  (fF-sec/kg) 

H  height  (m) 

Lm  mean  beam  length  (Eq.  22) 

M  molecular  wt 

M  average  molecular  weight  of  species  in  the  fire 
lit  mass  flow  rate  (kg/sec) 

ritr  fuel  evaporation  rate  (kg/sec) 

q  power  (W) 
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()grn  fllhySH c 

s  parameter  (Eq.  39) 

T  temperature  ("K) 

V  volume  (m3) 

Xo2  oxygen  mass  fraction  in  the  air 
r  radius  (m) 

Z  modified  Froude  number 

Greek  Symbols 

a  entrainment  coefficient 

j8  efficiency  of  combustion 

y  stoichiometric  oxygen/fuel  ratio 

f  emissivity 

8  defined  in  Eqs.  (10),  (19)  and  (30) 

k  absorption  coefficient  (m'~1) 

3  parameter  (Eq.  36) 

H  convergence  criterion  (Eq.  35) 

p  density  (kg/m3) 

a  Stefan-Boltzman  constant  ( W/m2  °K) 

r  transmissivity 

<l»  Fourier  coefficient  (Eq.  38) 

0  frequency 

w  defined  in  (Eq.  6) 

Subscripts 

a  environmental  air 

h  fuel 

c  combustion 

/  combustion  zone 

p  plume  zone 

r  radiation 

s  fuel  surface 

ss  combustion  zone — plume  interface 

t  flame  tip 

v  vapor 

O2  oxygen 

fuel  fuel 

0  flame-fuel  interface 

oo  ambient  conditions 

Superscripts 
0  base  value 
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